Non— Zero Fermi Level Density of States for a Disordered d-Wave Superconductor in 

Two Dimensions 



u 

< 



> 

o 

On 



I 

o 
o 



K. Ziegler*^, M.H. Hettler^ and P.J. Hirschfcld'' 
°' Max- Planck- Institut fiir Physik Komplexer Systeme, Auflenstelle Stuttgart, Postfach 800665, D-70506 Stuttgart, Germany 

Dept. of Physics, Univ. of Florida, Gainesville, FL 32611, USA. 
(February 1, 2008) 

It has been known for some time that, in three dimensions, arbitrarily weak disorder in uncon- 
ventional superconductors with line nodes gives rise to a nonzero residual density of zero- energy 
quasiparticle states N{fi), leading to characteristic low-temperature thermodynamic properties sim- 
ilar to those observed in cuprate and heavy- fermion systems. In a strictly two-dimensional model 
possibly appropriate for the cuprates, it has been argued that A'^(O) vanishes, however. We perform 
exact calculations for d- and extended s-wave superconductors with Lorentzian disorder and similar 
models with continuous disorder distribution, and show that in these cases the residual density of 
states is nonzero even in two dimensions. We discuss the reasons for this discrepancy, and the 
implications of our result for the cuprates. 

PACS numbers; 74.25-q, 74.25.Bt, 74.62.Dh 



Introduction. A good deal of evidence |0] has accumu- 
lated recently suggesting that the order parameter in 
the cuprate superconductors vanishes linearly at lines on 
the Fermi surface. Frequently these experiments have 
been interpreted in terms of a d,j.2_y2 pairing state, but 
states of "extended-s" symmetry with nodes have also 
been considered. One of the most interesting conse- 
quences of such nodes in three spatial dimensions is the 
creation of a nonzero density of zero-energy quasiparti- 
cle states A''(0) for infinitesimal disorder. Such a 
residual density of states (DOS) is of course reflected 
in many experimental observables, and may be shown 
Q to lead in particular to a term in the London 
penetration depth A(r) — A(0), and more generally to 
low-temperature thermodynamic and transport proper- 
ties characteristic of a normal Fermi system with strongly 
reduced DOS. Systematic .Z'n-doping and electron dam- 
age experiments have been found to lead to precisely 
these types of temperature dependences in YBCO sin- 
gle crystals [||, and in certain cases quantitative fits ||] 
have been obtained to the "dirty d-wave" model, in which 
the effects of potential scatterers on a 2D d-wave super- 
conductor are calculated using a t-matrix approximation 
assuming large electronic phase shifts. 

Recently, Nersesyan et al. Q have questioned the accu- 
racy of the t-matrix approximation when applied to a 
strictly 2D disordered d-wave system, pointing out that 
in 2D logarithmic divergences in multisite scattering pro- 
cesses, some of which are neglected in the t-matrix ap- 
proach, prevent a well-controlled expansion in impurity 
concentration. These authors avoid perturbation theory 
by using bosonization together with the replica trick, and 
predicted a power law DOS N{E) ~ a ~ 1/7, for 

sufficiently small energy E and disorder, rather than the 
analytic behavior N{u!) ^ const -f aE^ expected in 3D. 
They also argued that a non-zero DOS at = 0, a quan- 
tity indicating spontaneous symmetry breaking, may not 
occur because of the Mermin- Wagner theorem [nO| . 



Although the physical systems in question are in real- 
ity highly anisotropic 3D systems, the possibility of a 
2D-3D crossover at low temperatures could conceivably 
invalidate some of the results of the usual "dirty d-wave" 
approach. This would render the description of the low- 
temperature transport properties of the cuprate super- 
conductors considerably more complicated even were the 
order parameter to correspond to the very simple 2D 
dx2_y2 form usually assumed. 

It is therefore of considerable importance to check the 
results of Ref. |^ by other methods. In this paper we 
show that for certain types of disorder, exact results can 
be obtained for the DOS of strictly 2D disordered su- 
perconductors. We show that for any disorder diagonal 
in position and particle-hole space, the DOS of a classic 
isotropic s-wave superconductor has a rigorous thresh- 
old at the (unrenormalized) gap edge A, as expected 
from Anderson's theorem ||ll|]. Within the same gen- 
eral method, we show that the residual DOS N{0) of 
a superconductor with line nodes (e.g. d- or extended 
s-wave) is nonzero for arbitrarily small disorder, in dis- 
agreement with Ref. We believe that the DOS in a 
disordered system is not an order parameter which be- 
longs to the class of order parameters covered by the 
Mermin- Wagner theorem. This is supported by the fact 
that a non-zero DOS occurs also in other tight-binding 
models (e.g., model for two-dimensional Anderson local- 
ization [|l^), which are described by a field theory with 
continuous symmetry. 

As exact results are only obtainable for Lorentzian disor- 
der, we discuss ways [ p^ of obtaining information on the 
effects of other distributions, including models where the 
randomness has a compact domain. Finally, we compare 
our results to those arising from alternative methods, and 
comment on possible origins of the current disagreement. 
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Density of states. Here, we introduce a general method 
of calculating exactly the DOS of a superconductor for 
certain types of disorder, motivated by the analysis of 
Dirac fermions in 2D. jl^ The BCS Hamiltonian is given 

by 



(1) 



which describes quasiparticles in the presence of the spin 
singlet order parameter A. The Ui are the Pauli matrices 
in particle-hole space. The disorder is modeled by taking 
fj, = fj,x £ls a random variable distributed according to a 
probability distribution P{^x)- 

The kinetic energy operator — is taken to act as 
V2*(a;) = ^'(a;+2ei)+^'(a;-2ei)-t-5'(x-f2e2)-t-*(a;-2e2) 
on a function "^{x) of the sites a; of a 2D square lat- 
tice spanned by the unit vectors ei and 62. Note this 
function involves displacements of two lattice sites rather 
than one, as would be the case in the simplest tight- 
binding representation of the lattice kinetic energy. For 
a system of fermions in the thermodynamic limit, the 
bare kinetic energy will then have a band representation 
quite similar to the usual tight-binding form, with no 
particular distinguishing features near the Fermi level. 
The reason for this choice will become clear below. It 
obeys, of course, the same global continuous symme- 
tries discussed for the model in Ref. The bilo- 
cal lattice operator A = A^^x' is taken to act as a c- 
number in the isotropic s-wave case, A^'(a;) = A^(a;), 
whereas to study extended pairing we define A<i\If(a;) — 
Ad [-^{x + ei) + *(a; - ei) ± ^'(a; + 62) ± *(x - 62)]. 

We consider the single-particle Matsubara Green func- 
tion defined as G{iE) = {iEao — H)~^. We are 
primarily interested in calculating the DOS N{E) = 
— ilm^g(Gii(A:, ii? E + ie)) where (...) denotes 
the disorder average. The problem now is how to per- 
form this disorder average over the probability measure 
P{fix)dfix of the random variable fj,x- Exact results for 
the disorder-averaged propagator in noninteracting sys- 
tems can frequently be obtained for Lorentzian disorder, 
P{/j,x)dHx = {'y/Tr)[{nx - A«o)^ + l'^]~^diJ,x, by exploit- 
ing the simple pole structure of P{fJ.x) in the complex 
fix plane, /ig is the chemical potential of the averaged 
system. For convenience, we set fiQ = 0. The averaged 
Green function is {G{iE)) = J JJ^ dfixP{fJ-x)G{iE; fix), 
which may then be trivially evaluated if G can be shown 
to be analytic in either the upper or lower half-plane. 

In a superconductor, the Green function depends on the 
random variable fXx via fix i iE, as a consequence of the 
particle-hole structure. Therefore, the averaging of G 
with respect to Lorentzian disorder is not trivially pos- 
sible. However, we will show below that it is possible 
to reformulate the problem so that G consists of terms 
which are analytic in one of the half planes. This allows 
us to perform the averaging of the Green function for 



Lorentzian disorder. 

Isotropic s-wave superconductor. We first assume a ho- 
mogeneous s-wave order parameter, neglecting the re- 
sponse of the superconducting condensate to the random 
potential. The Matsubara Green function may be writ- 
ten G{iE) = -{iEao + H){E^ + H^)'^ where we note 
that — (— — fi)^<Jo + A^CTo since in the isotropic 
s-wave case, (— — /i)cr3 anticommutes with Acti even 
for random fi due to the locality of the order parameter. 

The expression + E"^ is proportional to the unit ma- 
trix; as a consequence, the Green function can be written 
in the simple form 



(-V2-^-iVA2+i?2)-l 

-(-V2-^ + ^v/A2 + ^2)-ll (2) 



It is straightforward to show that the imaginary part of 
this expression (after analytic continuation, iE — s- E+ie) 
for any given configuration of impurities is vanishing for 
\E\ < A. Therefore, the DOS shows a gap of size A 
independent of the distribution function P{fi). Thus, 
our model reproduces the famous Anderson theorem |Tl| ] 
which states that the thermodynamics of an isotropic s- 
wave superconductor are not affected by diagonal, non- 
magnetic disorder. The situation is different if the order 
parameter itself is random . 

d- and extended-s symmetry superconductors. The sec- 
ond class of examples includes the d-wave and extended- 
s "bond" order parameters A^ defined above. The cor- 
responding pure systems in momentum space fulfill the 
condition A^ = 0, so that nonmagnetic disorder must 
cause significant pair breaking ||]. The behavior of the 
imaginary part of the Green's function can be studied 
using a method analogous to that used for the s-wave 
case. However, the main difference is that the nonlocal 
order parameter term A^ai does not anticommute with 
(— — /i)(T3 anymore if fi is random. 

This requires a different type of transformation. We in- 
troduce a diagonal matrix (or staggered field) Dx^x' — 
[—\Y^+^'^5x,x' (note D'^ is the unit matrix). Now we 
may write 

= HDajDH ^ [(-V^ - fi)Dc7Q - iA^Dcr2] 

[D{-V^ - fi)(7o + iDAda2] (3) 

Because D commutes with — (as defined above) and fi, 
but anticommutes with the order parameter A , we have 
simply = H^, with H = (-V^ - fi)Dao - iAdDa2. 
Therefore, the quantity H"^ + E'^crt^ = {H + iEcro){H - 
iEao) can be used to write 



G{tE) = lilE^^ ( - iEaoy^ 



-{H + iEao) 



(4) 
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Observe that both H and H appear in this expression, 
but H only in the numerator. Since {H ± iE(To)~^ is an 
analytic function of /xD in either the upper or lower fiD- 
half plane and P{^D) = P{n) we can now perforin the 
disorder integration. The disorder averaged Matsubara 
Green function is translational invariant. Performing a 
spatial Fourier transform wc replace — by ^ = e^: — /io 
and find 



{G{^E)) 



{iE + ij)aa + ^0-3 + A ^ di 
(£; + 7)2+^2 + (A^)2 



GiiE^ 



(5) 



This is the Matsubara Green function of the pure system 
with the frequency iE shifted by the disorder parame- 
ter, iE ^ iE + ij. It should be noted that for the local 
(isotropic) s-wave order parameter discussed before the 
average over a Lorcntzian distribution in Eq. |2| implies a 
shift iVA^ + iV^^+E^ + ij. 

To obtain the DOS for the d-wave case we approxi- 
mate the sum over the momenta k in standard fashion 
as No '^^ lo^ ^ where No is the density of states of 
the normal metal at the Fermi level, with the tetragonal 
Fermi surface approximated by a circle. The result is 



N{E) = A^o 



2ir 



^^Im 
27r 



E 



(A2(0)-(£; + Z7)2)l/2 



(6) 



where the 

by Arf(0) 
27 



d-wave order parameter is approximated 
= AdCos(2(?i). At £; = 0, iV(0) = 
iVo;^ln(4Ad/7) for 7 « A^. Thus, the DOS is 
nonzero at the Fermi level for arbitrarily small disorder. 
For small E, N{E) varies as E'^. 

For more general continuous distributions P{^)d^i the 
averaged DOS can be estimated using again the ana- 
lytic structure of G. Applying the ideas of Ref. [p^ , 
one can derive a lower bound by a decomposition of 
the lattice into finite sub-squares. The average DOS on 
an isolated sub-squares can be estimated easily. More- 
over, the contribution of the connection between the sub- 
squares to the average DOS can also be estimated. A 
combination of both contributions leads to {N{0)) > 
ci min_^j<^<^j P{fi)^ where ci and fii are distribution 
dependent positive constants. In particular, fii must 
be chosen such that the spectrum of H{iiq = 0) = 
— V^CTs — iKdDai is inside the interval [— For 
all unbounded distributions, like the Gaussian distribu- 
tion used in Ref. , as well as compact distributions with 
sufficiently large support this estimate leads to a nonzero 
DOS at the Fermi level. 

Discussion and comparison to other methods. The ma- 
jor result in the d-wave (extended s-wave) case with 
Lorentzian disorder is the presence of a finite purely 
imaginary self energy Eg = — «7cro due to nonmagnetic 
disorder which leads to a nonzero DOS at the Fermi level. 



The latter is in qualitative agreement with standard theo- 
ries based on the self-consistent t-matrix approximation 
[ Apl as well as with exact diagonalization studies in 2D 
[16||. In contrast to such theories our self energy has no 
dependence on A^, i.e. it is the same as in the normal 
state. In Fig. 1 we show a comparison of the self energies 
of our theory and the limits of the t-matrix approxima- 
tion. 
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FIG. 1. Imaginary part of the self energy vs. frequency. 
For Lorentzian disorder (solid line) the self energy is constant 
17. The self energy of the self-consistent t-matrix approx- 
imation in the unitary scattering limit (dashed-dotted line) 
behaves cx (5A)^''^ at zero frequency. For Born scattering 
(dashed line) the value at zero frequency is nonzero, but ex- 
ponentially small. We have adjusted the impurity concentra- 
tion to obtain equal normal state self energies for the t-matrix 
results. 

A drawback of the model with Lorentzian disorder is that 
impurity concentration docs not appear explicitly in the 
theory. Whereas in the t-matrix approach we have with 
the impurity concentration and the scattering strength 
(or phase shift) two parameters associated with disorder, 
in the present model we have only 7, the width of the 
Lorentzian. A way of making a connection is by com- 
paring the variance of the Lorentzian distribution (7) 
and the variance of the distribution underlying the t- 
matrix approximation, which is a bimodal distribution 
of a chemical potential /i = /io with probability \ — 5 
(5 being the dimensionless impurity concentration) and 
M = Mo + ^ with probability 6 {V being the scattering 
potential). The variance Var^ of this distribution is de- 
termined by 



Var^ = (M^)-(/i)2 = y2(<5-0 



(7) 



For small concentrations of impurities, 5 « 1, we find 



Var. 



1/(5^/^. The 6^/"^ behavior is also found for 
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IniEo(i? = 0) in the t-matrix approach for strong scatter- 
ing. Since in our model the variance of the distribution 
is also the imaginary part of the self energy, this suggests 
that our model is closer to the strong scattering limit of 
the t-matrix approximation than the Born limit. 



discussions with R. Fehrenbacher and C. Mudry. They 
are also grateful for the support of the Institute for Fun- 
damental Theory, University of Florida, and the Institut 
fiir Theorie der Kondensierten Materie, Universitat Karl- 
sruhe. 



Finally, we comment on the discrepancies between our 
result and the calculation of Nersesyan et al., who found a 
power law for the averaged DOS with Gaussian disorder. 
One might question the analysis of Nersesyan et al. be- 
cause of the use of the replica trick, which is a dangerous 
procedure in a number of models, [ p^f However, Mudry et 
al. |l8| have obtained identical results for the continuum 
problem of Dirac fermions in the presence of a random 
gauge field using supersymmetry methods. We therefore 
believe that the crucial difference between our results and 
those of Ref. occurs in the passage to the continuum 
and concomitant mapping of the site disorder in the orig- 
inal problem onto the random gauge field. Only in the 
continuum case is there a direct analogy between disorder 
in the chemical potential and a gauge field; on the lattice, 
gauge fields and chemical potential terms enter quite dif- 
ferently. First, chemical potential terms are local while 
gauge fields are defined on bonds. Furthermore, chemi- 
cal potential disorder enters linearly in the Hamiltonian 
while gauge fields enter through the Peierls prescription 
as a phase in the exponential multiplying the kinetic en- 
ergy. 

Disorder of the gauge field type is furthermore nongeneric 
even in the continuum, as discussed by Mudry et al., 
who showed that the critical points of the system with 
random gauge field are unstable with respect to small 
perturbations by other types of disorder. pO| We expect 
that a proper mapping of the lattice Dirac fermion or d- 
wave superconductor problems to continuum models will 
inevitably generate disorder other than random gauge 
fields. Therefore, we believe that our result of a finite 
DOS at the Fermi level is the generic case for a d-wave 
superconductor in two dimensions. 

Conclusions. We have computed the single particle 
Green function and DOS for a model of a superconductor 
with nonmagnetic impurities. For an isotropic s-wave su- 
perconductor, we recover standard results; in particular, 
Anderson's theorem is reproduced. Our calculations for 
the disorder-averaged d- and extended s-wave propaga- 
tors show the DOS is nonzero for all energies, provided 
the distribution of the chemical potential is continuous 
and of sufficient width. The disorder average has been 
performed exactly in the case of a Lorentzian distribu- 
tion. Our calculation casts doubt on the result by Ners- 
esyan et al., who found a power law for the averaged DOS 
with Gaussian disorder, and suggests that the standard 
t-matrix approach to disordered d-wave superconductors 
is qualitatively sufficient. 
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It should be mentioned that even in the case of a random 
gauge field it is not established ]l^ , |l9t that the DOS in- 
deed vanishes at the Fermi level. 
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